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ABSTRACT 
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taken  to  obtain  the  correlation  distance  and  the  mean  square 
of  the  particle  velocity  using  the  autocorrelation  function. 
Finally  a  computer  program  was  written  to  evaluate  the  inte- 
gral wave  equation  by  numerical  methods.   It  was  found  that 
the  critical  parameters  in  the  mathematical  model  were  the 
correlation  distance  and  the  frequency  limits  of  integration. 
Small  variations  in  the  correlation  distance  modified  greatly 
the  width  of  the  predicted  beam  pattern,  while  changes  in  the 
limits  of  integration  had  a  moderate  effect.   The  Frequency 
Spectrum  was  obtained  in  the  anechoic  chamber  and  it  was  used 
to  determine  the  limits  of  integration  of  the  integral  solu- 
tion for  the  intensity  field. 

The  predicted  beam  patterns  were  in  satisfactory  agreement 
with  the  measured  far  field  beam  patterns.   The  differences 
between  them  are  probably  due  to  inaccurates  in  determining 
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I.   INTRODUCTION 

It  has  long  been  recognized  that  the  best  means  of  locat- 
ing a  deeply  submerged  submarine  is  by  detecting  the  sound 
that  it  radiates  into  the  water;  also  all  surface  ships  are 
subject  to  attack  by  passive  homing  torpedoes  and  acoustic 
mines,  both  of  which  operate  using  the  radiated  noise  of  the 
ship.   These  facts  make  the  far  field  radiated  signatures  a 
very  important  factor  to  determine  detectability  characteris- 
tics and  capacity  of  survival  in  a  hostile  environment. 

Free  field  measurement  in  the  far  field  involves  serious 
difficulties  like  ambient  noise,  multipath,   convergence  zone, 
attenuation,  among  others,  which  makes  impracticable  such 
kind  of  measurement.   Instead  a  prediction  of  the  Far  Field 
using  Near  Field  data  is  suggested.   Certain  general  princi- 
ples are  repeatedly  invoked  in  prediction.   The  first  is  that 
the  vibratory  powers  of  several  sources  is  additive.   This  is 
based  upon  the  assumption  that  the  signals  from  any  two 
sources  are  statistically  uncorrelated,  an  assumption  which 
is  almost  invariably  found  correct.   Signal  contributions  from 
non-overlapping  frequency  bands  are  inherently  uncorrelated; 
hence  the  total  power  over  all  frequencies  is  found  by  summing 
contributions  from  various  frequency  bands. 

Several  works  have  been  done  about  prediction  of  Far 
Field  beam  pattern  from  Near  Field  measurements  for  discrete 
frequencies  (Refs.  1,2,3,4).   In  general  these  methods  deal 


with  the  amplitude  and  phase  information  taken  at  the  Near 
Field  and  the  final  result  shows  that  the  Fourier  Transform 
of  those  measurements  lead  to  the  beam  pattern.   Consider 
the  random  noise  as  the  random  summation  of  sinusoidal  waves, 
thus  a  complete  analysis  of  the  frequency  bands  can  be 
achieved  using  a  very  narrow  filter  which  center  frequency 
can  be  shifted  to  cover  all  the  spectrum,  or  a  bank  of  narrow 
filters  can  be  used.   Adding  the  contribution  of  each  fre- 
quency it  should  be  possible  to  obtain  the  final  beam  pattern, 
but  either  procedure  involves  a  tedious  work,  specially  for 
broad  band  signals.   There  also  remains  the  question  of  cor- 
relation between  different  bands. 

The  more  general  problem  of  predicting  the  Far  Field  noise 
beam  pattern  from  a  random  noise  source  is  of  great  practical 
importance.   The  theory  of  the  noise  field  for  systems  that 
satisfy  the  wave  equation  has  been  developed  by  some  inves- 
tigators working  in  the  fields  of  acoustics  (Ref.  5),  and 
optics  (Ref.  6).   Different  methods  had  been  suggested  in  or- 
der to  generate  a  prediction  in  the  far  field,  one  of  them 
consists  of  using  correlation  functions  and  the  Green's  func- 
tion for  the  wave  equation,  for  general  time  variations.  This 
technique  is  used  by  Morse  and  Ingard  (Ref.  7)  in  their  book 
"Theoretical  Acoustic,"  for  describing  non-periodic  vibra- 
tions or  random  processes. 

This  paper  discusses  one  possible  method  of  Near  Field 
analysis  using  the  spatial  correlation  function  and  the 
Green's  function  which  enables  differential  equations  and 


boundary  conditions  to  be  combined  into  an  integral  equation 
for  the  wave  motion;  the  solution  of  the  integral  equation 
provides  the  intensity  in  the  Far  Field.   This  method  could 
practically  be  used  in  collecting  submarine  radiated  acoustic 
signatures . 


II.   GENERAL  THEORY 

The  analysis  of  the  disturbed  region  has  to  be  done  us- 
ing a  different  approach  than  for  discrete  frequencies  and 
one  of  the  most  common  becomes  the  correlation  technique. 
It  is  well  known  what  the  relationship  is  between  the  auto- 
correlation function  of  random  signal  and  the  total  average 
power  and  the  D.C.  power  contained  in  the  signal.   In  a 
sense  the  autocorrelation  function  is  a  measure  of  time 
variation,  spatial  variation  and  statistical  dependence  in 
the  region  of  inhomogeneities . 

The  study  of  the  statistics  of  noise  gives  the  possi- 
bility of  finding  the  solution  of  the  problem  using  the  cor- 
relation function  instead  of  a  frequency  decomposition  and 
individual  analysis  of  the  frequency  components  in  the  ran- 
dom signal.   In  the  process  any  phase  information  is  lost, 
but  the  important  point  becomes  the  relation  between  the 
autocorrelation  function  and  the  field  intensity,  which  is 
phase  independent,  and  the  angular  dependence  with  respect 
to  the  acoustic  axis  of  the  source. 

It  is  helpful  to  begin  with  a  non  periodical  vibration 
problem,  with  several  degrees  of  freedom,  each  degree  of 
freedom  has  its  own  equation  of  motion  and  the  general  equa- 
tion becomes: 

Mn  — «£■  =  F  (r)  +  F*  (r,t)  +  D  (v,t)        (2.1) 

dt 

n  =  1,2,. ..N 
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where   F  (r)   is  the  internal  force,   F'  (r,t)   is  the  exter- 
n  n 

nal  force,   D  (v,t)   is  the  frictional  force,  and   M   is  the 

'    n   '  '       n 

effective  mass  of  the  system. 

Plotting  the  behavior  of  one  component  of   x   and  its 
associated  velocity  on  the  Nth  phase  plane,  it  can  be  shown 
that  the  contour  is  not  a  closed- loop.   Usually  it  will  be  a 
complicated  non  periodic  curve  confined  within  some  limits 
given  by  the  energy  of  the  system  but  never  repeating  itself. 
Of  course  the  N  gegree  of  freedom  implies  N  different  frequen- 
cies in  such  a  way  that  the  combination  is  not  periodic.   The 
displacement  component  will  not  be  a  periodic  function  of 
time,  so  the  Fourier  Analysis  can  not  be  applied.   The  average 
value  of  the  displacement  is  also  time  dependent  for  short 
periods  of  time,  but  if  the  time  is  long  enough  the  fluctua- 
tion can  be  evened  out,  and  the  resultant  average  becomes  a 
constant;  again  it  is  necessary  to  use  the  statistical  pro- 
perties of  the  autocorrelation  function,  which  is  defined  as 
follows : 

R(x)    =    Lim    [i  fhT   x(t)    •    x(t+x)    dt]  (2.2) 

i-35T 

The  value  of   R(0)   is  equal  to  the  average  power  of  the 
signal  over  the  interval,  and  this  will  be  the  maximum  value 
of   R(t)  .   The  size  of  "t"  is  very  important  since  the 
rapidity  with  which   R(x)   drops  to  zero  as   t   is  increased 
is  a  measure  of  the  degree  of  randomness  of  the  oscillation 
of  the  function.   In  the  case  of  the  spatial  autocorrelation 
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function,  the  correlation  distance,  (d) ,  follows  the  same 
behavior. 

The  second  step  will  be  to  study  the  radiation  character- 
istic from  a  plane  surface  with  randomly  vibrating  area,  the 
goal  in  this  section  is  to  try  to  obtain  a  Green's  function 
such  that  one  or  more  parts  of  the  resulting  integral  equa- 
tion can  be  eliminated,  and  the  remaining  part  easy  enough  to 
be  integrated.   The  solution  of  this  problem  will  be  approached 
in  four  steps : 

(1)  Calculation  of  the  far  field  due  to  a  point  source  in 

the  presence  of  a  rigid  boundary  serves  to  outline  the  geometry 
and  describe  the  role  of  the  boundary  and  facilitate  the  choice 
of  the  right  Green's  function  for  subsequent  steps. 

(2)  Calculation  of  the  far  field  due  to  harmonic  motion  of 
parts  of  the  surfaces  itself  which  is  assumed  to  have  surfaces 
characteristic  of  a  rigid  boundary. 

(3)  Calculation  of  the  far  field  due  to  harmonic  motion  of 
parts  of  the  surfaces  which  have  a  particular  acoustic 
impedance. 

(M-)  Application  of  results  from  (3)  to  randomly  vibrating  area. 

1.   Consider  sources  in  the  presence  of  an  infinite  bar- 
rier like  the  arrangement  shown  in  Figure  2.1.   If  the  plane 
is  perfectly  rigid,  the  boundary  condition  is  that  the  par- 
ticle velocity  and  the  gradient  of  the  pressure  at  the  surfaces 
is  zero,  thus,  the  effect  of  the  boundary  plane  can  be  replaced 
by  a  set  of  images  sources,  symmetrically  located  with  respect 
to  the  plane,  both  radiating  into  the  region  containing   Po  . 
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The  acoustic  pressure  at   Q  ,      defined  by  the 
vector  r  ,   for  a  simple-harmonic  source  is  then: 


,    „      /     1         ikri  1         ikr?,.    -iu>t  /0    _n 

p   =   -ipok  Su   tj^r  e  +  Tf^r  e        2)e  (2.3) 


where  k  is  the  wave  number  and  S  represent  the  strength 
of  the  source  expressed  in  the  frequency  domain.   Some 
approximations  have  to  be  made  in  Equation  (2.3),  assuming 
that  Q  is  at  a  distance  such  that  the  two  waves  combine 
to  form  a  source  that  looks  like  a  single  nonsimple  source 
located  at  the  origin.   The  magnitude'  of  vectors  r-,   and 
r~   is  defined  by: 


rl  =   l^*^0l   =  r  "  Q'^o     =  r  "  ro  cos  9 


and 


rn  =   r-r1   =  r  +  2z^  cos  a  -  r  cos  <j> 

2     '    o '  o         o 


Using  these  approximations  Equation  (2.3)  becomes: 


.  „   .  ik(r  -  r   cos  6  +  z^  cos  a) 

_  -ipck  0     /,         n  o  o 

2  — 7tz —  Sf,  cos(kz^  cos  a)  e 

2ttt   a)    v   o 


The  distance  from  the  source  to  the  plane  is  an 

important  parameter  since  it  determines  the  constructive  or 

destructive  interference.   For  example,  if  kz   and  kr  are 

'      o       o 
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much  less  than  1  ,   the  source  is  much  smaller  than  a 
wavelength  and  the  behavior  of  the  pressure  in  the  far  field 
is  very  similar  to  that  produced  by  a  simple  source  of 
double  strength  of  the  original  source.   On  the  other  hand, 
if  z   is  equal  to  a  quarter  wavelength  the  separation 
between  the  source  and  its  image  is  one  half  the  wavelength 
and  its  effect  is  a  destructive  interference.   The  same 
effect  can  be  produced  if  instead  of  a  point  source  a 
multipole  source  exists  for  which  each  point  has  its  own 
image.   Of  course  the  boundary  produces  distortion  in  the 
near  field  which  is  not  considered  in  this  paper. 

2.   Now,  considering  that  the  sound  waves  in  the  bounded 
medium  may  be  produced  by  the  motion  of  the  plane  surface 
itself.   With  the  plane  rigid  at   z   equal  to  zero,  some 
portion  of  the  plane  may  be  moving  in  the   z  direction. 
For  a  simple  harmonic  motion  u  =  u  exp(-iwt)  ,   the 

Z      CO 

boundary  condition  is  that  the  gradient  of  the  pressure  is 
ikpcu  (x,y)exp(-iojt)  ,   whereas,  for  the  rest  of  the  plane 
the  partial  derivative  of  the  pressure  with  respect  to  z 
is  zero. 

The  time  independent  part  of  the  solution  of  the 
wave  equation  in  terms  of  the  Green's  function  is  given  by: 


Pu(r)  =  ///fu(ro)Gu(^)dvo  +  ffZQiS. )*   p  (r  )  .  p  (p  )fc  q ji)  ]  « 

O  O    O  0      0 


(2.4) 
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since  there  are  no  sources  in  the  medium  the  volume  integral 
is  equal  to  zero  and  the  surface  integral  involves  p  and 
the  outward  pointing  normal  gradient  at   z  =  0  ,   but  the 
description  of  the  physical  situation  has  specified  only 
the  normal  gradient.   The  Green's  function  should  be  selected 
in  such  a  way  to  eliminate  the  second  part  of  the  surface 
integral,  keeping  only  the  specified  normal  gradient  of  the 
pressure.   The  Green's  function  satisfying  this  requirement 
is  given  by: 


n    ikrn     ,    ikr0 

G  C — ;   =  -n e     +  t e 

o         1  2 


in  this  case  r,  -  r2  ,   because  the  source  is  in  the  surface 
boundary,  thus,  the  bounded  medium  Green's  function   (G  ) 
is  twice  as  much  as  the  unbounded  medium  Green's  function 
g   ,   the  reason  is  that  each  moving  element  of  area  radiates 
only  into  the  half  space   z  >  0  ,   the  source  and  its  image 
have  coalesced  to  form  a  source  of  double  strength. 

Applying  this  concept,  Equation  (2.4)  becomes: 

PuCx,y,z>  =  -ikpc  //uu(x0,y0)  VV^V  dX°  dy°     (2<5) 
=  -2ikpc  //uu(x0>yo)  S^3E^77)  dxo  dyo 


3.   A  perfectly  rigid  boundary  is  unrealizable.   The 
common  situation  is  that  the  boundary  has  a  particular 
acoustical  impedance  defined  by  its  density  and  the  velocity 
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of  sound  in  the  medium.   Then,  the  boundary  reacts  to  an 
incident  pressure  and  every  portion  of  the  plane  is  moving, 
and  a  different  method  has  to  be  applied  in  order  to  dis- 
tinguish the  active  regions  in  the  plane,  where  energy  is 
fed  into  the  medium  from  the  passive  areas. 

For  the  passive  areas  the  relationship  between  p 
and  its  normal  derivative  is  given  by: 


-u 


■  l&  -  A  P   atz  =  0  (2.6) 


z       ikpc    dz  pc 

where   (3   is  the  impedance  ratio  between  the  impedance  of 
the  medium  where  sound  propagates  and  the  impedance  of  the 
surface  boundary.   For  the  Active  Areas  the  normal  fluid 
velocity,  plus  the  admittance  of  the  boundary  times  the 
pressure  at   z  =  0  must  differ  from  zero,  this  difference 
will  be  defined  as  the  driving  velocity   (u  ) 


Cu  +-fi-p)     =   (J-iE+lp)     =  u  e~ia,t   (2.7) 
y    z        pc  p;z=0     kikpc  3z    pc  p;z=Q      u 


To  simplify  the  surface  integral,  the  Green's  function 

could  be  arranged  such  that  the  surface  integral  over  the 

passive  areas  is  zero.   This  is  accomplished  by  a  Green's 
function  that  satisfies: 


5^  ■  -lkBGu  (2.8) 
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The   integral  over  the   active   area  is 


pu     "     "  "CH-  Go>  "  PU  Ti*>      n     dxo^o  «.9) 

o  o    z=0 


Using  Eq.    (2.8)    in  Eq.    (2.9) 

P<o    "     "  "4lT  +  lk*Pj  Gw  dxo  dy0  (2-10) 

o 

The  expression  in  parentheses  is  the  driving  velocity  in 
the  active  region  given  by  (2.7),  thus  Eq.  (2.10)  becomes: 


pu  =  -ikpo  //uu(x0,yo)  Gu(^i^)  dx  dy0       (2.11) 

0   0 


which  has  the  same  form  as  Equation  (2.5),  but  this  time  an 
exact  representation  of  the  Green's  function  as  was  stated 
previously  is  not  possible.   However,  a  good  approximation 
valid  when  the  measurement  point  Q  is  not  closer  than  a 
half  wavelength  to  the  boundary  surface  and  which  satisfy 
Equation  (2.8)  is  given  by: 

ikr-,  ikr2 

G      .  t +    [    COS    a!    I    1    ]•    S (2.12) 

a)        Hirr-,  cos   a'    +    3  47rr2 

where  the  variables  correspond  to  Fig.  2.1.   Although  this 
form  of  the  Green's  function  is  derived  for  the  problem  of 
a  volume  source  in  the  presence  of  a  non-rigid  boundary,  it 
may  be  used  since  it  accomplishes  the  simplification  of  the 
problem  in  hand.   [In  this  approximation  there  is  an  image 
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at  P*  ,   but  the  wave  it  radiates  has  an  amplitude  which 
depends  on  the  angle   a'   which  the  line   rp  makes  with 
the  z-axis  normal  to  the  boundary.   The  ray  of  acoustic 
radiation  which  reaches  the  point  Q  from  the  source  point 
at  P  by  reflection,  strikes  the  reflecting  surface  at   Q 
at  an  angle  af   of  incidence  and  reflection.   The  wave 
which  goes  directly  from  P   to   Q  without  reflection, 
given  by  the  first  term,  is  independent  of  the  inclination 
of  r-.  ;   the  reflected  wave,  given  by  the  second  term, 
depends  on  the  angle  of  reflection.   The  factor  3  introduced 
the  variation  in  the  reflected  sound  due  to  the  boundary.   In 
this  way  the  term  in  parentheses  represents  the  reflection 
coefficient. ] 

In  the  far  field  Eq.  (2.12)  can  be  simplified  using 

rl  =  r2: 


ikr1  ikrn 

g       1  -to  1 

G 


e r_    cos  a'  -  B-i  _  e r   cos  a'   -.     ,p  ,-x 

4Trr1   L1   cos  a!  +  3 J    27rrx   Lcos  a'  +  SJ    ^d'1^>> 


From  Eqs.  (2.11)  and  (2.13)  the  pressure  in  the  far 
field  is  given  by: 

ikr. 


Observe  that  the  integral  is  over  the  active  region 
of  the  plane  which  is  assumed  to  be  confined  to  the  circle 
of  radius  "a"  about  the  origin. 
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In  the  far  field  the  last  equation  can  be  simplified 
neglecting  terms  of  order  a/r  and  smaller  and  using  the 
approximation   cos  a'  =  cos  a  =  z/r  , 

r.,   ~   r  -  (x  cos  <J>  +  y  sin  <p)    sin  a 


—   =  —      and    tan  d>   z  - 
r,     r  Y   x 


Thus,  Equation  (2.7)  becomes: 


ikr  -ik*W  sin  a 

p  (r)  z        -ikpc  r  *ros  a —  •  | ffu   (x  ,y  )e     °     dx  dy 

*w  w      8  +  cos  a   2-rrr     oo   o,,7o  o  J o 

(2.15) 


where  the  vector  w  lies  in  the  xy  plane,  with  x  ,  y 
as  component;  and  vector  k  has  the  r  direction  whose 
components  are: 


k   =   k  cos  <f)  sin  a 

A 

k   =  k  sin  ({>  sin  a 

k   =   k  cos  a 
z 


A  close  study  of  the  integral  part  of  Equation  (2.15) 
shows  that  the  angular  distribution  of  the  Far  Field  is 
proportional  to  the  two  dimensional  Fourier  Transform  of  the 
distribution  of  the  driving  velocity  in  the  active  region, 
which  is  defined  by: 
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WV      =     Z?       ""        VVo'    e"V0"V°     dxo  <^o 


ik  x  +k  y 
VW     =     "         WV   e  dkx  dky 

—  00 


2 
fu>(a>(,))      =      gVcos^x      Vk   cos    ♦    sln   a>    k   sin   ♦    sin   a)       (2-16) 


Thus  Equation  (2.15)  takes  the  form: 

ikr 
Pu(r)   S   -ikpc  ^r-     tJoL,*)  (2.17) 

The  angular  dependence  is  given  by  the  factor   f   ,   which 

CO 

measures  the  relative  amplitude  of  the  pressure  with  respect 

to  the  center  of  the  active  region. 

Several  properties  of  the  Par  Field  can  be  studied. 

If  the  coordinate  axes  are  rotated  to  X  .Y  ,  such  that   X 

o*  o'  o 

axis  is  the  intersection  of  the  rz  plane,  the  appropriate 
transformation  is  given  by: 


xo  =  Xo  cos  *  "  Yo  sln  * 


yo  =  Xo  sin  *  +  Yo  sin  * 


Under  this  transformation  the  time  independent  particle 
velocity  becomes: 


21 


,  -ikX    sin  a 

U  (k  cos  <t>  sin  a,k  sin  <J>  sin  a)  =  -~  f  e  dX    /  u  (X  ,Y  )  dY 

0)  liTT  WOO      O 


This  equation  shows  the  one  dimensional  Fourier 
transformation  in  the  X  direction,  of  the  average: 


u.(o),X  J   =  i-  /  u  (X  .Y J  dY^ 

4>'0        2tT      0)0   0      o 


of  the  particle  velocity  u  (x0jy0)   in  the  Y   direction. 

It  can  also  be  observed  that  u.(o),X  )   is  a  function 

of  X   ,   the  distance  in  the  direction  of  the  measured 
o 

point  Q  ;   it  is  also  a  function  of  $  ,   the  direction  of 
the  rz  plane,  and  of  the  frequency   o)/2tt  . 

If  the  particle  velocity  u  (x,y,t)   is  an  arbitrary 
function  of  space  and  time  and  u  (x,y)   is  its  time  dimension 

0) 

Fourier  transform,  then  the  pressure  field  at  a  point  Q  (r,a»<j>) 
is  the  inverse  transform  of  P   given  by  Equation  (2.17) 

-io)(t  -  |) 
p(r,t)   =  ^   /   f/cctj))  e  do)        (2.18) 


Now,  applying  the  rotation  of  axes  as  defined  before, 
Equation  (2.18)  can  be  expressed  as: 

-lo)(t  -  f)  °°    i(|)XQ  sin  a 

p(r,t)     =    *g-  •  s  ?os  a  /    -lo)e  do)     /  e  dX 

^     *  2irr      g  +  cos  a  o 

—CO  -oo 

CO  00 

x    /     dY      i     /    e^r    u(X  .Y  ,Z)  dZ 

o     2tt  oo  (2>19) 


p(r,t)     =    £-  ■  g  ;os  a        |r    [/  dX   /dY    u(X.Y.t  +  i  X    sin  a  -  £)] 
•  2rrr      $  +  cos  a     9t  o        o        o'  o'         c     o  c 
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This  equation  shows  clearly  some  of  the  properties  of  the 
Far  Field.   An  appropriate  summation  of  the  normal  acceleration 
of  the  active  area  gives  the  total  field  at  the  point  of 
interest.   Also  the  contribution  to  the  pressure  at  r  ,   at 
time  t  ,   from  the  line  X  ,   -»  <  Y  <  »  ,   depends  on 
the  normal  acceleration  on  the  line  at  time 

t+(l/c)X  sin(a)-r/c  ,   at  time  r/c  -  (X  /c)  sin(a)   earlier, 
which  is  the  time  required  for  the  sound  to  travel  from  the 
line  to  point   Q  . 

4.   All  the  formulae  developed  so  far  can  be  used  to 
compute  the  mean  intensity  at  r,a,<J>,  produced  by  an  area  A 
of  the  boundary  which  is  vibrating  with  a  random  distribution 
of  velocity.   The  particle  velocity  in  the  active  area  A  is 
a  function  of  position  as  well  'as  of  time.   In  this  regard, 
it  becomes  a  random  variable  and  a  way  to  describe  its 
behavior  is  using  its  statistical  properties.   In  this 
particular  case  the  autocorrelation  function  is  used  as  a 
means  to  describe  the  behavior  of  the  random  source. 
Considering  the  particle  velocity  as  a  function  of  position 
and  time  in  the  active  area,  the  autocorrelation  function 
is  defined  by: 


RU.ti.t)   =  ^r  ///  u*(x,y,t)  uz(x+^,y+n,t+x)  dx  dy  dt 


where  the  integration  is  over  the  active  area  A  ,   and  T 
is  over  a  long  interval  of  time.   Also  the  Fourier  transform 
of  the  velocity  is  given  by: 
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,  iajt-ik   x-ik  y 

U    (k    ,k    )      =     -±*  fffu   (x,y,t)    e  x  y      dx  dy  dt    . 

oj      x     y  Q^5  z 


From  this  equation  the  space  time  spectrum  level  of 
the  particle  velocity  is: 


p  at  iuT-ik j-ik  n 

U  I*     ■       -^TT    fff     R(£>n>T)e  y  d^dndx  (2.20) 


In  the  far  field  it  is  known  that  the  intensity  is 
given  by  the  real  part  of  the  pressure  times  the  particle 
velocity,  but  also  the  velocity  is  related  to  pressure  for 
the  relation  u  =  p/pc  ,   thus  using  Eq.  (2.16)  and  Eq.  (2.17) 

«        u    2,  2    2  2 

l-E-r  a  ^£  ^-^ — — — %   |U  (k  sin  a  cos  <j>,  k  sin  a  sin  <j>)  I 
M       r    | 3  +  cos  a| 

(2.21) 

This  equation  is  related  to  the  autocorrelation 
function  by  Equation  (2.20).   The  problem  is  now  reduced  to 
getting  the  form  of  R(5,t\»t)   in  space  and  time.   A 
reasonable  assumption  could  be  a  Gaussian  waveform  where  the 
value  of  R  falls  off  smoothly  as  its  arguments  are  increased 
Assume  that  the  autocorrelation  is  given  by: 

1,-2.2*    1  ,    n2 
RU,n,T)   =   |UQ|2  e  2d  (2f.22) 


2 
where   U    is  the  mean  square  particle  velocity  of  the  active 

area  averaged  over  time  interval  T  ;   d   stands  for  the 

correlation  distance  along  the  plane,   oo   is  the  upper 
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angular  frequency  limit  of  the  spectrum,  which  also  can 
give  the  correlation  time,  and  the  lower  limit  of  the 
wavelength. 

Using  this  value  of  the  autocorrelation  function 
and  equations  (2.20)  and  (2.21),  the  intensity  field  was 
evaluated: 

flT           ,/2    UV|«  |2cos2  a   -¥f  >2  ^  «  "  kr>2     <2-23> 
T   /    n   ^         AT  /0  \-5/2              '    o1  _  n 

Ir(w)  =  pc  —  (2tt)   ^        -p ^  e 

r         c  a)  1 3  +  cos  a  I 

n'       ' 

But  this  is  only  for  a  specific  frequency;  in  order 
to  get  the  total  field  intensity  it  would  be  necessary  to 
consider  the  contribution  of  all  the  frequency  components 
of  the  random  signal.   A  logical  way  to  accomplish  this  is 
to  integrate  Equation  (2.23)  with  respect  to  frequency. 
Note  that  the  integrand  has  a  maximum  when  the  exponent 
is  zero  at: 


a)  =  21/2  un(l  +  kn2d2  sin2  a)1/2 


A   .  lUol2cos2a      (knd>2 


r       _«,  2irr    iB  +  cosap   [1  +  (k^draiiT  a]^ 

(2.24) 

A  simple  analysis  of  this  formula  shows  that  the 
intensity  is  proportional  to  the  impedance  of  the  medium, 


25 


times  the  mean  square  particle  velocity  in  the  active  area, 
times  the  ratio  between  the  active  area,   A  ,   and  the  area 
of  the  half  sphere  of  radius  r  ,   times  an  angle  distribution 
function  ,  which  depends  on  the  maximum  frequency  of  the 
spectrum  (inverse  of  the  correlation  time),  the  speed  of 
sound  and  the  correlation  distance,   d  . 

Trying  to  apply  this  theory  to  practical  cases,  two 
main  problems  arose.   The  first  is  that  in  real  situations 
the  frequency  spectrum  is  band  limited,  usually  due  to  the 
physical  characteristics  of  the  source.   Under  this 
consideration  the  limit  of  integration  should  have  an  upper 
bound  (UL)  and  a  lower  limit  (LL),  in  this  form,  Equation 
(2.23)  becomes: 

TOO  TTT  1/UEL2      .      2     .       1/    U)n2      , 

tm    2    M  rar**    d 1U°'  cos  a      2*   /      2-*-)slna-2^  * 

I  (r)  =  2pc  -p  (2tt)   -"         -s P   •  -m-     /       u    e 

r  ceo  16  +  cos  al  TT 

n1  ■  LiL 

(2.25) 

There  is  no  analytical  solution  for  this  equation,  so 
a  numerical  method  of  integration  had  to  be  used  in  order  to 
get  the  intensity  in  the  Far  Field. 

The  second  problem  is  the  difficulty  in  measuring  the 
autocorrelation  function  of  the  randomly  vibrating  area, 
which  is  related  to  the  particle  velocity  in  the  active  area. 
Instead  it  can  be  assumed  that  if  the  area  close  to  the 
active  region  is  vibrating  with  similar  characteristics, 
then  the  impedance  ratio  becomes  unity  and  the  autocorrelation 
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in  the  Near  Field  gives  the  mean  square  particle  velocity; 

Equation  (2.25)  becomes: 

1,0x1%  2  .2  1,  uk2   , 

-,2|TT    i2  2  UL  2    c  2  ui 

A            -V2     d   'Uol      C0S     a  2  n 

I  (r)  =  2pc  -\  (2ir)  ^    -g— 2 ?      /    w2e 

r  c  a:  1 1  +  cos  a  |       TT 

n  ^  (2.26) 


This  equation  is  angle  dependent,  thus  the  intensity 
in  the  Far  Field  depends  on  the  angular  position  with  respect 
to  the  acoustical  axis.   In  this  way  it  is  possible  to 
evaluate  the  beam  pattern.   Observe  also  that  u)   should 
lay  inside  of  the  limit  of  integration  or  at  least  equal  to 
the  upper  limit,  otherwise  the  process  would  show  a 
discontinuity  which  is  not  possible  because  the  frequencies 
contributing  to  the  field  are  given  by  the  limit  of 
integration. 
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III.   EXPERIMENTAL  RESULTS  USING  THE  PLANE  MODEL 

To  validate  the  theoretical  model  described  in  Section 
II,  Far  Field  and  Near  Field  measurements  were  made  on  a 
horn  speaker  in  an  anechoic  chamber.   The  actual  beam  pattern 
in  the  far  field  of  the  horn-speaker  was  taken  and  also  the 
autocorrelation  function  in  different  points  of  the  near 
field  were  obtained  in  order  to  get  the  mean  square  particle 
velocity  in  the  Near-Field  and  be  able  to  predict  the  Far 
Field  Beam  Pattern. 

A.   FAR- FIELD  MEASUREMENT 

The  normalized  Far  Field  beam  pattern  was  obtained  using 
the  experimental  set  up  shown  in  Figure  3.1.   Since  the 
source  was  driven  by  random  noise  the  distance  of  the  near 
field  is  variable.   The  criterion  for  near  field  is  usually 
given  for  one  particular  frequency,  as  the  frequency  is  made 
lower  and  lower  the  near  field  also  decreases  and  eventually 
disappears  for  practical  purposes.   That  is  the  reason  for 
the  presence  of  the  high  pass  filter  which  limits  the  low 
frequency  component  to  1300  Hz  in  order  to  get  an  appreciable 
near  field  region.   This  is  necessary  for  obtaining  far  field 
and  near  field  measurements  under  the  same  frequency  spectrum, 
The  signal  is  then  amplified  and  applied  to  the  horn  speaker. 
The  microphone  was  placed  in  the  far  field  at  a  distance  of 
508  meters  from  the  source,  the  received  signal  was  amplified 
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and  the  output  connected  to  a  true  RMS  Voltmeter  which  also 
generates  a  DC  voltage  output  proportional  to  the  deflection 
of  the  needle  for  the  selected  scale,  this  voltage  was  ap- 
plied to  the   Y   axis  of  the  plotter.   The   X   axis  of  the 
plotter  received  information  from  the  position  indicator  as 
a  DC  voltage  proportional  to  the  angular  displacement  of  the 
acoustic  axis  of  the  horn  speaker.   The  normalized  pressure 
amplitude  versus  angular  position  is  shown  in  Figure  3.2. 

This  was  the  actual  Far  Field  Beam  Pattern  which  had  to 
be  compared  with  the  beam  pattern  evaluated  using  the  mathe- 
matical model  described  in  the  previous  section. 

B.   NEAR-FIELD  MEASUREMENT 

The  intensity  in  the  Far  Field,  as  defined  by  Equation 
(2.26)  depends,  among  other  factors,  upon  the  mean  square 
particle  velocity.   The  mean  square  particle  velocity  is  a 
source  quantity  which  was  shown  to  be  related  to  the  auto- 
correlation function.   The  evaluation  of  the  factor  is 
possible  by  using  an  array  of  microphones  over  the  entire  ac- 
tive region  and  the  summation  of  the  received  signal  fed 
into  the  correlator.   The  microphone  array  was  not  available, 
instead  a  single  microphone  was  used  to  take  a  measurement 
each  5  cm  in  the  vertical  and  horizontal  direction  trying  to 
cover  the  active  area  as  shown  in  Figure  3.3.   The  autocorre- 
lation values  obtained  at  each  single  point  is  then  averaged 
and  the  resultant  value  was  used  to  describe  the  mean  square 
particle  velocity. 
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The  correlator  was  a  SAICOR  model  SAI-4  3A,  which  was 
used  in  the  autocorrelation  mode,  the  function  is  determined 
simultaneously  at  400  incremental  lag  points,  the  completed 
correlation  function  is  displayed  at  one  time,  the  function 
is  provided  with  incremental  lag  or  time  delay  value  rang- 
ing from  0.2  microseconds  to  1  sec  resulting  in  total  time 
delays  from  80  microseconds  to  400  sec.   Precomputation  de- 
lay of  800  lag  values  selected  in  200  lag  increments  allows 
the  autocorrelation  function  to  be  viewed  symmetrically 
about  zero  or  up  to  800  lag  values  removed  from  zero.   The 

averaging  is  accomplished  digitally  in  a  linear  fashion  with 

9      17  . 
fixed  summations  ranging  from  2   to  2    in  binary  steps. 

The  quantization  is  8  bits  per  channel. 

The  SAI-43  operates  as  follows:   one  input  signal  is 
applied  to  either  channel,  and  internally  routed  to  the 
second  channel,  the  signal  is  sampled  at  the  rate  indicated 
by  the  operator-selected  sample  increment  and  then  is  quan- 
tized in  both  channels.   Channel  I  is  the  channel  which  con- 
tains the  400  word  memory.   Channel  II  contains  the  one  word 
memory  -  the  present  sample.   At  the  time  called  "Present 
Time" ,  the  correlator  performs  the  multiplication  of  the 
present  value  word  of  Channel  II  by  each  of  the  400  words 
in  Channel  I.   The  product  "0"  times  "0"  becomes  an  entry  in 
R(0),  the  first  point  of  the  autocorrelation  function  the 
product  "0"  times  "1"  becomes  an  entry  for  R(l)  and  finally 
"0"  times  "399"  becomes  an  entry  for  RC399).   At  the  end  of 
this  arithmetic  cycle  the  correlator  has  performed  one  calcu- 
lation for  each  so-called  bin  of  the  autocorrelation  function 
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The  product  in  R(0)  is  the  first  entry  for  the  zero-lag  auto- 
correlation value,  while  the  product  in  R(39  9)  is  the  first 
entry  for  the  correlation  value  at  a  lag  in  Channel  I  equal 
to  399t.   These  first  entries,  resulting  from  one  sample  in 
Channel  II  are  expanded  to  the  full  correlation  function  by 
taking  repeatedly  updated  value  of  the  sample  in  Channel  II 
and  correspondingly  updated  values  in  Channel  I  and  then 
averaged.   In  effect  the  400  word  memory  slides  ahead  one 
word  and  the  arithmetic  cycle  begins  again.   For  all  t  incre- 
ments that  are  equal  or  longer  than  5  00  usee  the  SAI-43A  is 
said  to  operate  in  real  time.   But  for  increments  less  than 
500  usee  the  figure  is  different.   Assume  an  increment  less 
than  50  msec,  it  could  have  completed  eight  cycles  in  400 
msec  instead  of  one,  therefore  it  is  eight  times  slower  than 
the  real  time  instrument. 

The  output  of  the  correlator  can  be  displayed  in  an 
oscilloscope  or  in  an  x-y  plotter,  the  voltage  output  in  the 
autocorrelation  mode  can  be  determined  by  the  following 
formula: 

R(v2)  =  2£gOLT).(10)dBA/10  (3>1) 

where   x   is  the  voltage  as  seen  as  the  output  of  the  SAI-43A 

dBA  is  the  dB  setting  of  the  input  channel.   R   is  the 

2 

correlation  value  in  (volt)   to  be  determined. 

For  illustration  purposes  Figure  3.4  and  Figure  3.5  show 
the  autocorrelation  function,  R( 0,0,0),  for  the  point  of 
maximum  intensity,  (center  of  the  speaker),  and  the  point  of 
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minimum  intensity,  which  correspond  to  the  corner  of  the 
boundary  of  the  active  region,  the  relative  power  between 
them  is  more  than  2  0  dB. 

The  physical  arrangement  in  the  anechoic  chamber  and  the 
instrumentation  for  taking  the  autocorrelation  function  at 
each  point  in  the  near  field  is  shown  in  Figure  3.6.   The 
autocorrelation  of  the  band  limited  noise  becomes  the  sin  x/ 
x,  function,  the  amplitude  of  the  function  is  variable  depend- 
ing on  the  position  of  the  microphone.   As  described  before 
it  is  also  possible  to  get  the  power  at  each  point  of  the 
autocorrelation  function. 

The  criteria  used  for  determining  the  active  area  was  on 
the  basis  of  the  differences  in  power  between  the  center  of 
the  horn  up  to  the  point  which  had  2  0  dB  less  power,  which 
corresponds  to  the  boundary  of  the  active  area  shown  in 
Figure  3.3. 

The  other  important  measurement  was  the  frequency  spectrum 
of  the  received  signal  which  is  given  in  Figure  3.7.   This 
was  obtained  using  the  basic  setup  but  with  a  wave  analyzer 
instead  of  the  correlator.   The  frequency  sweep  was  made  very 
slow  in  order  to  get  a  good  average.   The  main  use  of  this 
figure  was  to  provide  the  proper  limits  of  integration  for 
Equation  (2.26).   The  limit  of  integration  became  a  quite 
important  parameter  since  it  modifies  greatly  the  form  of  the 
normalized  beam  pattern;  thinking  in  terms  of  a  transducer, 
it  is  known  that  higher  frequencies  make  the  beam  narrower 
and  the  low  frequencies  have  a  wider  beam,  so  looking  at 
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Figure  3.2,  the  actual  beam  pattern,  it  is  logical  to  assume 
that  the  main  contributions  are  given  by  the  low  frequencies, 
and  that  is  exactly  what  the  Frequency  Spectrum  shows  in 
Figure  3.7.   The  lower  limit  of  integration  is  assumed  to  be 
1100  Hz  and  the  upper  limit  is  5000  Hz,  the  contributions  of 
the  other  frequencies  was  small;  this  also  showed  that  the 
speaker  has  a  low  frequency  response,  although  it  covers  the 
important  band  of  the  human  hearing. 

Another  critical  factor  in  the  Near-Field  measurements 
was  the  correlation  distance,  since  its  effect  in  the  width 
of  the  beam  is  great.   The  symmetry  of  the  source  caused  some 
difficulties;  it  was  found  that  all  points  located  at  the 
same  distance  from  the  axis  of  the  horn  had  a  very  good  cor- 
relation, regardless  of  the  distance  between  them,  but  along 
the  same  angular  position  a  distance  of  12  cm  was  enough  to 
diminish  greatly  the  spatial  autocorrelation  function.   This 
suggested  that  the  correlation  distance  for  this  experiment 
is  a  vectorial  parameter,  but  was  used  like  a  scalar  along 
this  work.   It  is  useful  to  point  out  that  in  a  real  situa- 
tion with  an  asymmetric  source  driven  by  several  sources  with 
different  frequency  bands,  like  a  submarine  in  which  each 
section  is  affected  for  different  set  of  frequencies  and  also 
the  physical  configuration  changing  in  space,  the  autocorre- 
lation distance  has  to  be  easy  to  measure  experimentally  and 
its  vectorial  characteristic  should  disappear  leading  to  a 
scalar  parameter.   Otherwise  the  pattern  must  be  built  up 
from  several  sources. 
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Along  this  experiment  a  correlation  distance  of  12  cm 
was  used. 

C.   RESULTS  FROM  NEAR-FIELD  MEASUREMENT 

The  first  step  in  the  data  processing  was  to  get  the 
average  mean  square  particle  velocity.   Using  the  data  col- 
lected in  the  experiment  and  the  formula  given  by  Equation 
(3.1)  it  was  possible  to  evaluate  the  power  at  each  point  of 
the  active  area,  this  power  was  averaged,  the  entire  process 
was  done  using  an  IBM-3  60  computer  program  in  fortran  IV 
language.   The  value  obtained  in  this  program  represents  the 
mean  square  of  the  particle  velocity  expressed  by  the 
Equation  (2.22),  where  the  value  of  the  autocorrelation 
function  is  at  zero  spatial  and  time  shift. 

The  approximation  used  for  that  evaluation  was,  of  course, 
an  additional  source  of  inaccuracy  in  the  experiment,  but 
should  not  affect  too  much  the  final  result,  because  what  is 
really  affected  is  the  intensity  field.   The  normalized  beam 
pattern  is  practically  unaffected  since  the  relative  power 
between  each  angular  position  was  the  subject  of  interest. 

In  the  case  of  trying  to  obtain  a  radiation  pattern,  which 

2 
is  beyond  the  scope  of  this  work,  then  Uo   becomes  a  very 

important  factor  since  it  affects  the  Power  Level  of  the 

signal. 

The  solution  of  Equation  (2.26)  is  the  actual  prediction 

of  the  Far  Field  Beam  Pattern.   Again  the  computer  system 

was  used  for  performing  numerical  integration,  a  library 

sub-routine  called  DQATR  was  used.   It  was  necessary  to  use 
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double  precision  in  order  to  obtain  a  reasonable  accuracy. 
The  sub-routine  uses  the  trapezoidal  rule  in  connection  with 
Rombergs  principles  for  evaluating  the  approximation  of  the 
integral  value  of  the  function. 

The  resulting  computed  Far-Field  is  shown  in  Figure  3.8 
in  comparison  with  the  actual  Far-Field.   Note  that  the 
agreement  between  the  predicted  and  the  measured  beam  pattern 
is  fair  in  shape  but  slightly  out  in  values,  the  actual  beam 
pattern  shows  a  wider  beam  that  the  predicted  one.   Several 
factors  can  produce  that  difference,  the  correlation  distance 
could  be  one  of  the  main  factors  due  to  its  specially  vari- 
able characteristic  in  this  experiment,  as  it  was  explained 
before.   Also  the  limits  of  integration  assumed  an  important 
role,  strong  high  frequency  components  could  produce  that 
narrow  beam. 

In  order  to  study  the  variation  of  the  function  with 
respect  to  the  limits  of  integration  more  calculations  were 
made  reducing  the  upper  limit  to  42  50  Hz,  all  other  parameters 
were  kept  constant.   The  result  of  this  calculation  is  shown 
in  Figure  3.9.   This  time  a  better  agreement  was  obtained 
between  the  actual  and  the  predicted  beam  patterns.   A  wider 
beam  was  obtained  but  still  showed  a  small  difference  with 
the  actual  beam  pattern. 

The  effect  of  the  correlation  distance  is  shown  in  Figure 
3.10.   The  limits  of  integration  were  kept  at  42  5  0  Hz  and 
1100  Hz,  and  the  correlation  distance  was  increased  to  15  cm 
in  the  calculation.   It  can  be  seen  that  a  small  variation  in 
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this  parameter  affects  greatly  the  predicted  width  of  the 
beam.   This  Figure  also  shows  that  the  beam  pattern  is  more 
sensitive  to  the  correlation  distance  than  to  the  limits  of 
integration. 
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IV.   CONCLUSION 

The  possibility  of  predicting  the  far  field  beam  pattern 
from  randomly  vibrating  surfaces  has  been  demonstrated.   It 
is  feasible  to  obtain  the  intensity  field  and  the  radiation 
pattern  from  near  field  measurements  of  distributed  sources. 
Although  the  theoretical  model  developed  in  Section  II  was 
tested  experimentally  with  satisfactory  results,  more  work 
needs  to  be  done  in  both  areas .   The  amount  by  which  the  beam 
pattern  depends  on  each  parameter  makes  it  difficult  to  state 
a  solid  argument  for  describing  the  range  of  validity  of  each 
one. 

The  computed  beam  pattern  was  obtained  solving  equation 
(2.26)  using  a  computer  program  for  solving  the  integral  by 
numerical  methods.   The  limits  of  integration  were  chosen  by 
means  of  the  Frequency  Spectrum,  which  was  taken  using  a  wave 
analyzer  with  automatic  frequency  sweep.   The  maximum  range 
considered  was  from  1100  Hz  to  5000  Hz.   It  was  assumed  that 
the  frequencies  producing  the  field  were  inside  of  this  fre- 
quency band. 

The  mean  square  of  the  particle  velocity,  which  has 
energy  information,  was  shown  to  be  related  to  the  autocorre- 
lation function.   A  correlator  was  used  for  that  measurement 
and  the  power  contained  in  the  signal  was  easily  calculated 
with  equation  (3.1).   Note  that  in  this  case  a  correlator  was 
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used  but  a  true  RMS  voltmeter  can  also  be  used  to  measure 
the  power  of  the  signal. 

The  correlation  distance  was  also  measured  with  the 
correlator,  a  pair  of  microphones  were  located  in  the  active 
area  and  the  two  dimensional  (spatial  and  time)  correlation 
function  was  observed.   A  correlation  distance  of  12  cm  was 
used  for  calculation,  but  the  measured  value  was  found  to 
vary  with  location  of  the  microphones.   The  calculated  beam 
pattern  is  very  sensitive  to  this  parameter.   Its  effects 
are  bigger  than  those  due  to  the  limits  of  integration. 

Different  configuration  of  sources  should  be  tested  ex- 
perimentally.  Sources  should  be  excited  for  different  fre- 
quency bands  in  such  a  way  as  to  produce  an  active  region 
with  a  well  defined  correlation  distance.   This  will  solve 
one  of  the  main  problems  of  this  experiment.   Also  it  is 
suggested  that  the  frequency  spectrum  of  the  configuration 
should  not  necessarily  be  continuous,  as  in  this  form  the 
arrangement  will  be  closer  to  real  cases. 

The  theoretical  model  does  not  consider  the  spectrum 
level  of  the  signal  for  the  different  frequencies.   This 
should  be  tested  more  carefully  in  order  to  state  a  rule  for 
the  limits  of  integration.   In  the  case  of  discontinuities 
in  the  Frequency  Spectrum  the  model  fails  to  give  a  straight- 
forward solution,  although  it  looks  reasonable  to  perform  a 
summation  of  the  contribution  of  each  frequency  band.   It 
was  not  possible  to  verify  this  concept  with  the  experimental 
setup. 
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Agreement  was  generally  good  between  measured  and  com- 
puted patterns,  especially  when  the  upper  limit  of  integra- 
tion was  reduced  to  425  0  Hz.   The  actual  beam  pattern  was 
slightly  broader  than  the  computed  pattern  in  all  cases. 
This  is  not  surprising  since  the  far  field  beam  pattern  de- 
pends on  the  frequency  response  of  the  source  and  the  limits 
of  integration  were  chosen  in  a  rather  empirical  way,  so  the 
computed  beam  pattern  could  include  a  bigger  frequency  band 
than  the  real  one  and  the  width  of  the  beam  would  be  too 
narrow. 

The  technique  developed  here  might  also  be  useful  in 
measuring  the  radiated  noise  pattern  of  submarines.   Some 
additional  considerations  should  be  taken  into  account,  and 
the  theoretical  model  should  be  validated  by  actual  measure- 
ment in  anechoic  tanks. 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  EVALUATING  THE  FAR  FIELD  BEAM  PATTERN 
USING  THE  THEORETICAL  MODEL  DEVELOPED  IN  SECTION  II 


C       PROGRAM  TO  EVALUATE  THE  BEAM  PATTERN  OF  A  HORN 

C      SPEAKER  DRIVED  WITH  RANDOM  NOISE. 

C       FOR  PROGRAMING  PURPOSES  THE  FOLLOWING  SYMBOLS  HAVE  BEEN 

C      USED: 

C       CI  =  CHARACTERISTIC  IMPEDANCE  OF  THE  MEDIUM 

C       D  =  CORRELATION  DISTANCE. 

C      BETA  =  ANGLE  WITH  RESPECT  TO  THE  ACOUSTICAL  AXIS. 

C      T  =  AVERAGE  OF  THE  MEAN  SQUARE  OF  THE  PARTICLE 

C      VELOCITY  IN  THE  ACTIVE  REGION. 

C       XL  =  LOWER  LIMIT  OF  INTEGRATION. 

C      WN=UPPER  LIMIT  OF  INTEGRATION. 

C       A  =  AREA  OF  THE  ACTIVE  REGION. 

C       C  =  SPEED  OF  SOUND  IN  THE  AIR. 

C 

C 

C       CALCULATE  THE  INTENSITY  IN  THE  FAR  FIELD 

C 

C 

IMPLICIT  REAL*8  (F-H,Q-Z) 

REAL  LABEL/'     V 

REAL*8  TITLE  (12) 

REAL*8  PI,A,AUX(999),D,C,CI,CS 

DIMENSION  P(73) ,BETA(73) 

EXTERNAL  F 

COMMON  D,C,RBETA,SN,PI 

PI  =  3.1^15926 

D  =  0.12 

EPS  =  0.001 

WN  =  2.*PI*5000. 

NDIM  =  999 

A  =0.5**2 

T  =  0.062 

CI  =  415. 

C  =  3^3. 

R  =  5.8 

XL  =  2.*PI  *1100. 

51  =  CI*A*  T  *(D**2) 
WP  =  4250.  *2.  *  PI 

52  =  (R**2)  *(C**2)  *WP*((2.*PI)**1.5) 
DO  10  I  *  1  ,37 

BETA  =  BETA(I)*PI/  180. 

SN  =  (DSIN(RBETA))**2 

CS  -  DCOS(RBETA) 

CALL  DQATR(XL,WN,EPS,NDIMSF,Y,IERSAUX) 

P(I)=(S1*(0S**2)*Y)  /  (S2*((1.+0S)**2)) 
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BETA  (74-1)  =  -BETA(I) 
P(74-I)  =  P(I) 

WRITE(6,100)  Y,IER,BETA(I),P(I) 
10  CONTINUE 
100  FORMAT(/'Y  = » ,El6 . 7 , 5X, « IER= ' ,15, 5X, 'BETA= ■ ,P7 . 2 , » INT* 

1  f,El6.7) 
C 
C 

C       CALCULATE  THE  NORMALIZED  AMPLITUDE  OF  THE  PRESSURE  FAR 
C       FIELD  BEAM  PATTERN. 
C 
C 

Z  =  P(37) 
DO  20  1=1,73 
P(I)=  ((P(I)/Z)**0.5) 
WRITE(6,200)  BETA(I),P(I) 
P(I)  =  P(D*  2. 
20  CONTINUE 
200  FORMAT(/, 'BETA  =  ' ,F7 . 2, 13X, 'NORMALIZED  AMPLITUDE  =  ',F 

17.5) 
READ(5,50)  TITLE 
50  FORMAT (6A8) 

CALL  DRAW (73, BETA , P , 0 , 0 , LABEL , TITLE ,20. ,0.5,0,4,0,2, 
19, 8,0, LAST) 
WRITE(6,400)  LAST 
400  FORMAT(/rLAST  =',I2) 
STOP 
END 


FUNCTION  F(W) 
IMPLICIT  REAL*8  (A-P,0-Z) 
COMMON  D,C,RBETA,SN,PI 
BIP  =0.5 
BFP  =  4250. 
ERP  =  2. 

WP  =  BFP  *  PI  *  BRP 

F  =(W**2)/((DEXP(BIP*((W*D/C)**2)*SN))*(DEXP(BIP* 
1((W/WP)**2)))) 
RETURN 
END 
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